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Theory for the Rabi and internal Josephson effects in an interacting Bose gas in the cold collision 
regime is presented. By using microscopic transport equation for the density matrix the problem is 
' mapped onto a problem of precession of two coupled classical spins. In the absence of an external 

, excitation field our results agree with the theory for the density induced frequency shifts in atomic 

■ clocks. In the presence of the external field, the internal Josephson effect takes place in a condensed 

{^jQ' Bose gas as well as in a non- condensed gas. The crossover from Rabi oscillations to the Josephson 

i~{ , oscillations as a function of interaction strength is studied in detail. 

^ ■ PACS numbers: 03.75.Fi, 05.30.Jp 
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^ ■ I. INTRODUCTION 

\ Recent experiments achieving the trapping and Bose condensation of atoms with hyperfine structure [^-Q have 
'J^' created a lot of theoretical interest in the behavior of atomic gases with internal degrees of freedom |]8|-|l5|] . Most of 
this interest has focused on how the presence of internal degrees of freedom effect the spatial behavior, such as the 
hydrodynamical modes or vortices. An equally important question posed by these experiments is how the internal 
states of the atoms are effected by interactions. 
^ , In a variety of experiments the transitions between internal states are used to either manipulate or probe cold 
gases. Fountain atomic clocks [p^-p^ use Ramsey spectroscopy technique to measure the energy differences between 
^ \ internal states very precisely, and have found that this energy difference depends on the densities of the different 
^ • internal states in the gas. Another class of experiments have used Electromagnetically Induced Transparency |2^] 
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techniques to excite the atoms into a superposition of internal states by absorbing a 'probe' beam and were able 
reemit the beam after some delay. These experiments measured the change in the 'stopped' light beam as the gas 
evolved in the superposition state . Both of these classes of experiments measure the internal dynamics of the 

atoms in a normal gas. Experiments involving change in the internal state of the atoms have also been carried out 
' ' ■ on Bose Einstein Condensates (BEC). In the MIT spin polarized Hydrogen experiment ||2^,Q, BEC was detected 
by changing the internal states of the atoms from IS to 2S using a two photon excitation. Two photon transitions 
have also been successfully employed by the JILA group in a number of experiments to create topological excitations 
p5| , p6t or phase textures p7| | in a Rubidium BEC. One can expect many related experiments to be carried out in near 
OO future as the study of cold atomic gases is a rapidly expanding field. 

Although all these experiments try to elucidate the internal dynamics of a cold gas, it is useful to classify them into 
^""^ , two groups with respect to the strength of the external field that couples to the internal states. In the weak field limit, 
' experiments probe the internal dynamics of a gas of atoms in the absence of an external field. For example, in the 
, atomic clock experiments, although strong fields at the beginning and the end of Ramsey spectroscopy are required, 
an essentially free evolution of the wavefunction in a superposition of two internal states is measured. In the MIT 



^ hydrogen experiment ||23|, the number of excited atoms is a small fraction of the total number, and thus all the effects 



measured are linear in the excitation rate. Thus the external field is only a weak probe of the free dynamics. On the 
other hand, in the strong field limit the internal dynamics is sensitive to the intensity of the external field, and the 
number of atoms changing their internal state form a large fraction that may reach 100%. In the JILA experiments, 
. ^ supplying a TT pulse all the atoms are transferred from one internal state to the other, and duration of the pulse is 

k> , the inverse of Rabi frequency, showing that the internal dynamics depends on the strength of the external field. It 
; I • should be also kept in mind that in all these experiments, even in the weak field limit, the external field has large 
' enough number of quanta to be treated as a classical field. Also we will assume throughout that the excitation is 
carried out by a coherent field. 

The interaction between two atoms generally depends on their internal states |^,^. For the dilute cold gases 
considered in this work, interactions can be simply characterized by the s-wave scattering lengths a^p which depend 
on the internal states a and (3 of the interacting atoms. Thus a two state problem will have three different scattering 
lengths. To be in the purely s-wave scattering limit the gas has to be cold enough, so that 
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where a is the largest of the s-wave scattering lengths in the problem and At = /i(27rmfc;,T)~^/^ is the thermal de 
Broglie wavelength. The temperature range defined by the criterion (|l|) corresponds to the so-called cold collision 
regime. 

It is important to realize that the inequality (0) not only defines the condition for the s-wave scattering to be 
dominant, but also is a criterion for quantum statistical effects (indistinguishibility of the particles) to be observable. 
To see this, consider an evolution of the state of one particular atom. This atom will experience two different kinds 
of scattering processes: the coherent forward and backward scattering which preserves correlation of the momentum 
state with the internal state, as well as ordinary elastic collisions causing scattering at an angle different from or tt 
which destroy this correlation |30|j3l| ]. The forward and backward processes occur with the 'rate' 

r-1 = ^n, (2) 



while the ordinary elastic collisions rate is 



^coH = ■i'^a^nvT, (3) 



where vt — \P^^k^Tjm is the thermal velocity. Comparing these two rates one can see that the coherent processes 
happen at a much higher rate in the temperature range (nl). It is then possible to measure the frequency shifts of 
the internal states and decide the statistics of the system. The gases satisfying the condition (|^) are called quantum 
gases pO|-|3^. 

For the densities relevant for the alkali gas experiments, the temperatures below which the gas becomes quantum 
in the sense (|l|) are much larger than the Bose-Einstein condensation temperature. [This is true because quantum 
degeneracy becomes important only at the lower temperatures for which the thermal de Broglie wavelength is the 
order of average interparticle spacing. At — n~^l'^ \ For Bose systems the substantial overlap between wavefunctions 
of different particles and bosonic tendency to be in the same state then causes Bose-Einstein condensation, where 
a finite fraction of the particles of the gas share the same spatial wavefunction. The internal states of the particles 
define the interactions between the particles, thus determine the self-consistent spatial wavefunction that all the 
particles share j^J^, transition between internal states are also effected in return. Most importantly, there are no 
backward(exchange) scatterings between two bosons sharing the same spatial wavefunction, and the mean field energy 
shift for a particle in the condensate is different from a particle in the normal gas. 

An important example of the modification of the transitions between the internal states upon condensation is 
given by the internal Josephson effect [p^- ^8|j29|] . In this effect, the internal states act as the two reservoirs and 
the external field provides the weak link between these two reservoirs similar to the tunneling coupling through a 
barrier in the ordinary spatial Josephson effect p9|-^. Although the JILA experiments have come close, the 
internal Josephson effect has not been observed in the alkali gases to date, mostly due to the fact that the shape of 
the condensate is different for the two internal states. Any transition between the two internal states causes spatial 
disturbances, making the internal Josephson oscillations hard to observe. It should however be possible observe the 
effect by either using a very shallow trap or going to the adiabatic limit in which the Rabi frequency for transitions 
is much larger than the oscillation frequencies of the condensate in the trap |^,^. We will consider the case of a 
homogeneous gas and investigate the Internal Josephson effect, having in mind a gas sample in a very shallow trap. 

To understand these different phenomena in a general framework, we first derive the transport equation for the 
density matrix of the system, treating the external degrees of freedom semi-classically. The transport equation can be 
used to understand transitions between internal states of a Bose gas in the presence of spatial non-uniformity caused 
by a trapping potential, as well as in the presence of macroscopic flows, such as currents created by a vortex. 

Below we focus on the dynamics of internal states in a spatially uniform system. We assume that all the excitation 
fields are applied uniformly on a homogeneous sample. We find the equations of motion for the internal state density 
matrix of both the normal component and the condensed part of a Bose gas. Then we use a Bloch sphere representation 
p6| for the density matrix to get an equivalent equation of motion for classical spins. 

After deriving the equations of motion we consider two simple cases: a normal gas and a fully condensed gas. 
We find the free precession frequencies for the relative phase of two internal states, which leads to an expression for 
the density dependence of the frequency shift in fountain atomic clocks. Our results for the normal gas agree with 
previous work p3|-^, and we present a new result for the fully condensed case. 

Then we investigate the response of the Bose gas to an external field when it is either above condensation temperature 
or at zero temperature, i.e. fully condensed. For the fully condensed Bose gas we show how the Rabi oscillations in 
a non-interacting gas turn into the internal Josephson oscillations, and calculate the oscillation frequencies exactly 
for the entire parameter range. By repeating the same analysis for a non-condensed gas, we find that the internal 
Josephson effect exists in a cold non-condensed gas. This somewhat surprising result shows that the Josephson effect 
does not require a broken symmetry, and can be observed in any coherent system. We calculate the frequencies for 
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the normal case and compare with the results for condensed gas. We also compare the internal Josephson oscillations 
with the spatial Josephson oscillations discussed in [^9|-^ and show that 'macroscopic quantum self trapping' effect 
is also present for the internal Josephson effect. 

In the next section, we study the internal dynamics of a Bose gas which has comparable densities of condensed and 
normal components. We first discuss the dynamics without an external coupling field and calculate the oscillation 
frequencies for the condensate density in each internal state. In the limit of small oscillation amplitude out result 
agrees with the resonance frequencies calculated for the MIT hydrogen spectrum |Q . We briefly discuss the behavior 
of the system under a strong external field. 

Finally we give the transport equation and the equations of motion for the internal density matrix of a Fermi gas. 
We consider a freely precessing state and its response to an external field, and compare the results with those for the 
Bose gas. 



II. THE TRANSPORT EQUATION 

A dilute Bose gas in a confining potential U{r) can be described by the Hamiltonian 



n = /dV^^+(r) ( + hU4r) ) Mr) (4) 



+ Id'rJ2^i^+ ir)^P+ {r)i^p (r)^„ (r) , 

a, 13 

where the Greek indices a, (3 run over all internal states. The potential energy hUa{r) includes the internal energies 
hwa of the states \a), hUa{r) = hU{r) + hwa , and the interaction parameters Aq^ are related to the s-wave scattering 
lengths aai3 as 

Xap — aap- (5) 

m 

For a problem with N internal levels there are N{N + l)/2 interaction strengths Aq,/3 to be specified. The operators 
'ip{r) are the canonical second-quantized Bose operators, satisfying 

[V'a(r),^^(r')]=5a/3<5(r-r'). (6) 

Transitions between internal states caused by an excitation field, are described by another term to be added to the 
Hamiltonian (Q): 

Hexc = Y.j d^V^p (r, t)^t (r^p (r) ■ (7) 

a/3 

We assume harmonic time dependence 

Vai3{r, t) = Aai3{r) exp{ina0t) (8) 

where Aa/^ (r) is a function of coordinates describing spatial distribution of the excitation field [ pO| . The requirement 
of being hermitian gives Va/^it) = V^^{t). 

In the case of Bose condensation we can treat the operators as having non-zero expectation values. This is 
achieved by a substitution 

Mr) ^ Mr) + Mr), (9) 

where ipa is a c-number representing the condensate wavefunction. 

To describe the dynamics of the internal states of the normal gas, we need the populations of each level, as well as 
the coherences between any two states. Thus we define the full density matrix g{r,p) and the internal density matrix 
p{r), both being N x N hermitian matrices for N internal states: 

ecpir^p) = J d^r'i^'tir + ^^)Mr " ^))e^''-'' (10) 

/(Pp 
j^Qa0{r,p) (lAa W^A/jW). 
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The internal dynamics of the condensate can be described by specifying complex amplitudes of internal states, in 
total N complex numbers for an N component condensate. We find it useful, in order to keep track of relative phases 
of condensate internal states, to introduce an N x N matrix similar to the density matrix of normal states pajsir)- 
We will also need to keep track of the condensate flow described by JajS- These quantities are defined as 



(11) 



Here one notes that p is always rank one, thus the matrix does not contain any additional information but is introduced 
just for convenience. The diagonal elements of paa{i~) and Paa{r) give the thermal and condensate populations nta(r) 
and Ucair), respectively. Total density of atoms in state \a) is na{r) = nta{r) + nca{r). 

We start with deriving the equations of motion for g and p. Using the Heisenberg equations of motion for the field 
operators 



(12) 



and then taking the required expectation values we obtain ^ (r, p) and ^ (r) . The contribution of each term in the 
Hamiltonian to the evolution of g and p can be evaluated separately. To remind the reader how this calculation is 
carried out we present the contribution of potential energy to ^{t,p)- First we calculate 



V'a(r)=i / dV^C/^(r')V^(r')^^(r'),V'a(r) (13) 

7 

= i j d\'Y,U,{r')[^p+{r%'^^{r)]^l^+{r') 
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Upon complex conjugation we get 

V'^(r)=iC/„(r)V+(r). (14) 
Then we write the contribution of potential energy to ^ (r, p) as 

£0°'\r,p) = J d'r'e^P-^' ((V'+(r + r' /2)Mr - r' /2)) + {^+{r + r'/2)^^{r - r' /2))) (15) 

~ i J rfV'e^f'-' (Uo^ir) - Upir) + r' • v,( ^"(^) + )^ (^+(r + r72)V/3(r - r'/2)) 

= i{U^{r) - U0{r))g^0{r,p) + v^( ^"(^) + ) . Vpg^0{r,p). 

Going through this procedure for each of the terms in the Hamiltonian we get the equation of motion. Whenever 
we need to calculate the expected value of a 4 particle operator we do it using Wick's theorem, 



(16) 



In this averaging the first term is referred to as the direct (Hartree) term as it depends only on densities while the 
second is the exchange (Fock) term that involves coherences. 

We find the following equations describing the space and time variation of the two matrices: 



(17) 



+ i 



Bai3{r,p) 
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■ i ^7" Pq7 {r)gji3 (r, p) - >'0jPj°p {r)gaj{r, p) 

7 7 

7 

7 7 

"*E [^T"^''' ^)BiP^'^^P) - ^/57(^ 0£'a7('''P)] 

7 

{r,p)] 



and 



+ 1 



E ^7"P"7('')P7/3(^) - E ^PlPll3i'^)Pci{r) 
7 7 



(18) 



where p*°^{r) = p(r) + p(r). 

Let us summarize here the approximations used to derive the above equations and assess their vaUdity. The Hartree- 
Fock approximation is used to treat interactions, that is the same as keeping the lowest order term in the expansion 
in the gas parameter na^ . This is an excellent approximation for the dilute atomic gases where na^ is usually about 
10-6. 

The most important restriction on the above equations ( [r^ , [T8| ) is that they are valid for times not longer than the 
elastic scattering time 



Trail = ^T^a^nVT. 



(19) 



However, for a cold gas the characteristic frequency of the evolution described by (|l^,|l8|) is much higher than r^;;, 
and this approximation gets better as the temperature is decreased. The decoherence due to elastic collisions can be 
described heuristically by adding a damping term 



• (Damp) A 2 ( ^ s: \ 

0af3 = -47ra„^'yT£'a/3(l - da/3) 



(20) 



on the right hand side of Eq. (|T^) . 

All other effects, such as the difference of the trap potential for different internal states, or density oscillations 
created by a change in the internal state are adequately described by (|l7|,|l8|) within the Hartree-Fock approximation. 



III. EQUATION OF MOTION FOR THE INTERNAL DENSITY MATRIX 

Our main interest in this paper will be the consequences of interparticle interactions in the internal state dynamics, 
so we will simplify the above general equations by assuming that we are in a uniform system, with no normal or 
condensate flows, and further assume that the external fields are applied uniformly. In this case we will have no need 
for the full density matrix g{r,p), as internal dynamics is completely described p{r) which will be the same at all 
points r. 

Thus equations of motion reduce to 

P77' = ^{'^■y ~ ^7' +U^f — U^i)p^^' + i E XajiPia + Pfa)Paj' (21) 

a 

AcK-y' (pcK-y' ~t~ /^a^' ) ~t~ ^ ( (t) p^^y' ^'y'a(t^P'ya^ 
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where are defined as, 



+ U^- Uy)p^Y +?^Aa7P7aPQ7' (22) 

a 
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The dynamics described in the above equations are fairly comphcated if there are fields coupling more then two 
states, so we will concentrate on the simplest case, where only two internal states are coupled, and all the populations 
in the other states remain constant. In this case we will be concerned with the dynamics of two, two by two matrices. 
One for the normal gas and one for the condensate. Let us assume that states 1 and 2 are coupled. We can go to a 
Larmor basis with the frequency of the coupling field so that the elements of our density matrix (and similarly the 
condensate matrix) can be redefined. The diagonal elements of the matrices do not change while the off-diagonal 
elements change as 

p'i2 = Pi2 exp[-irJi2i], = (p'la)*- (24) 

In this basis we will find it useful to rewrite the equations of motion in the Bloch representation |46| . Being 2 by 2 
hermitian matrices we can expand, 

P^^, = PoS-yY + S ■ CT-y^' (25) 

where a are the Pauli matrices. 

In this representation po and po will be proportional to the total number of atoms in the condensate, and above the 
condensate respectively. The z components of both spins are proportional to the population difference between the 
two internal states, belonging to the normal gas or the condensate. For the normal part, the norm of the projection 
of spin onto the x — y plane represents the degree of coherence between the two internal states. For the condensate 
spin however, this norm is directly proportional to the geometric mean of the populations of the two internal states. 
The angle corresponding to a rotation around the z axis is related to the relative phase of the two internal states, in 
both cases. 

The matrix Vap has no diagonal elements as defined in the excitation Hamiltonian. However, we can take the 
detunings that come as a result of going to Larmor basis, as the diagonal elements. Vii(V22) being defined as 
+ (— )[w2 — Wi — fli2]- Then we can again expand V^^' as 

V^^' = VoS^^i + V ■ a^ij. (26) 
In this representation we can write the equations of motion for the condensate and the normal gas spins as : 

S ^ S X Bn + 2Xi2S X Sc + 2S xV (27) 
Sc = ScX Bc + 2Xi2Sc X S + 2ScXV 

with, 

Bn = [(All - A22)(2po + Po) + (All + A22 - 2Ai2)(25 + 5,) • z]z (28) 
Be = [(All - A22)(/5o +Po) + (All + A22 - 2Ai2)(S^+5c) • z]z. 

It is important to observe that the equations conserve the total densities in the condensed and non condensed 
fractions of the gas. To get a better sense of the equations we can note that they can be derived from the XXZ self 
interacting Hamiltonian: 



n = (All - A22)(po + Po)z ■ (S + Sc) + (All - A22)po5 • S (29) 
-2{S + S,).V 



I jij Qi oj _|_ jij oi oj I jij Qt oj 
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where, 



2Ai2 

J = I 2Ai2 I (30) 
(All + A22) 



using the Poisson spin algebra: 



{S\S^ = e'^'^S'' (31) 
{S\Si} = 0. 

Let us briefly remind what the physical bases for various terms in the Hamiltonian are. The factor of 2 difference in 
the normal spin-normal spin interaction compared to the condensate spin-condensate spin interaction stems from the 
fact that there are no exchange scatterings within the condensate. The terms coupling the condensate spin with the 
normal spin arise from processes in which an atom in one internal state goes through an exchange scattering process 
with a condensate atom in another internal state. 



IV. DENSITY SHIFTS IN ATOMIC CLOCKS: FREE PRECESSION OF A SINGLE SPIN 

We start the analysis of the dynamics described in Eq. (|27| ) by considering the situation in which we can represent 
the system by a single spin. There are two such cases. The first is at zero temperature, when almost all the atoms 
are in the condensate; while the second is above the transition temperature, when there is no condensate present. 
However, our analysis is still restricted to temperatures low enough to satisfy the quantum gas condition. At > as- 

We first analyze equations of motion Eq.(|2^) in the case where there is no external field. Such an analysis is needed 
to understand any interaction related effects in experiments where particles spend a substantial amount of time in a 
superposition state, such as atomic clocks p6|-|l9t. 

Let us first consider the case where there is no Bose condensation. Thus our equations ( |27| ) will be reduced to 

S = Sx Bn. (32) 

We can easily see that the z component of the spin will be conserved and we will only have a precession around the 
z axis and the precession frequency will be 

u) = (w2 - wi) + 2(A22n2 - Aiini + Ai2(ni - ^2)) + ^ (A2'y - Ai-y)n^. (33) 

This result agrees with the theory of frequency shifts in atomic clocks [^3|-^5| . 

We can now ask the same question for a sample which is fully condensed. When all the atoms are in the condensate, 
they will all be sharing the same spatial wavefunction, this will mean that no exchange processes will be possible. 
That would eliminate the factor of 2 multiplying the combination of A and ri's for states one and two in the previous 
equation. Not surprisingly we get, for a fully condensed sample 

UJ ^ {w2 ~ wi) + {\22n2 - Aii^i + Ai2(ni - 712)) + ^ (A2^ - Ai^)n^. (34) 

Although it is not desirable to use a condensate in a fountain atomic clock due to very high density shifts that 
would result. We see here that such an experiment can be used to probe the correlations in a condensate by simply 
looking at the density induced frequency shifts. 



V. INTERNAL JOSEPHSON EFFECT: SINGLE SPIN UNDER AN EXTERNAL FIELD 



In this section we consider the response to an external field, only when the system can be represented by a single 
spin, i.e. the system is either fully condensed or not condensed at all. In the former case, when the system is at zero 
temperature we have two condensates in internal states 1 and 2, which are connected by the "weak link" supplied by 
the external field. This is the well studied problem of the internal Josephson effect [|4[-Q . 
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The important difference between the spatial Josephson problem and the internal Josephson problem is that in 
the latter atoms in different internal states interact, while in the former two particles in different reservoirs do not. 
However if the interaction between the two internal states is taken to be zero 

Ai2 = 0, (35) 



then equations for the internal Josephson problem reduce to that of the spatial problem |}40|-[42[. thus the internal 
effect displays all the phenomena the spatial Josephson effect has. 

As in the spatial Josephson problem there are three different regimes depending on the strength of the coupling 
between the internal states [ p9| . If the external field is much stronger than the interactions 

\V\ > \n, (36) 

the system is said to be in the Rabi regime. (Here A can be taken to be the largest of X^p and n the total density.) In 
this regime interactions are not important and the system experiences large oscillations in the density of each internal 
state. In the opposite regime, when the external field is extremely weak, 

\V\ < \n, (37) 

the system is in the Fock regime. The amplitude of density oscillations is negligible and the phase difference between 
the two states evolves as in the free precession case discussed in the previous section. Between these two regimes is 
the Josephson regime, for which, 

\V\ - Xn. (38) 

In this regime, the density imbalance between the two internal states go through small oscillations, similar to the 
usual spatial Josephson effect. 

All three regimes are described by Eq.|^^, which reduces to 

Sc^Sc^Bc + 2Sc X V (39) 
Be = [(All - A22)po + (All + A22 - 2Ai2)5c • z]z, 

at zero temperature. 

As the dynamics conserves the magnitude of the spin we only have two dynamical variables, which can be taken as 
angles {9, (f) in the spherical polar coordinates representing the orientation of the spin. Another conserved quantity 
in this dynamics is the Hamiltonian introduced earlier, 

n = (All - A22)poi • + \j''SlSi + 2Se ■ V. (40) 

On the sphere defined by {6, cj)), contours of constant H will define the paths along which the spin will precess. We 
can write the Hamiltonian as a function of {9, (f) as 

n = ylcos2(6') + Bcos(6') +Csin(6')cos((/)) + D, (41) 

and identify 

^=^(Aii+A22-2Ai2)|5c|' (42) 

((All -A22)po|^c|+21/^|^e|) 

C = 21/^1^,1 
D^Xi2\Sc?. 

At strong external fields the paths are almost circular with centers on the line oriented along vector V , passing 
through the origin, as in the usual Rabi problem, and significant changes in the populations occur throughout the 
course of a Rabi oscillation. As a function of 9 and cf) the Hamiltonian has one maximum and one minimum, thus all 
the trajectories circle these extremum points, if the external field satisfies 



C_ 
2A 



>(1 



B 
2A 



2/3 



N3/2 



(43) 
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When the intrinsic Rabi frequency \V\ becomes the order of the density caused shifts ~ An, We get into the 
Josephson regime. Even on resonance, large population transfers from one state to another does not take place. This 
can be easily seen from the structure of trajectories on the Bloch sphere. When the external field does not satisfy 
(^), instead of having one minimum and one maximum, the Hamiltonian has two maxima, one minimum and a 
saddle point. The two constant H trajectories crossing at the saddle point separate the sphere into three regions, 
giving us three kind of trajectories, circling around the minimum or one of the two maxima. None of these oscillations 
however, result in large population transfers, which is in contrast with the Rabi problem, in which one can change 
internal states of all the atoms with an arbitrarily small field on resonance. Josephson effect can be understood if 
we realize that the transition frequency for an atom depends on the populations. So for weak fields even a small 
population transfer carries the transition away from resonance, i.e. makes the effective detuning much larger than 
the Rabi frequency. 

We can calculate the oscillation frequency along each path, for both high and low fields. The period of the precession 
along a path C on which H{9, (p) — Ti' is given by 

By converting the integral to a surface integral over a 6 function and integrating over the angle cj), we have 

or 1 
T = — / (45) 
\A\ J Va;4 + Csx^ + Cax^ + Cix + Co 



where 



Ca = 2^ (46) 
; C^ + B^ H_D 

B{H-D) 



Ci = -2 



C^-jH- Df 
l2 • 



In the high field case, for the external field satisfying (^), the Hamiltonian will take values between Hmax and 
Hmin, the values of the Hamiltonian at the maximum and minimum points. Any value of H in this range will uniquely 
correspond to one trajectory and the period of motion on such a trajectory will be given by 




^_ 4 1 {xi-X2r-ip-q)\ 

^-MV^^(2V M ^ ^''^ 

~ (m — .xi)^ 
= (m ~ X2f' 

where K is the complete elliptic integral p7[ |, x\, X2 are the real roots, m is the real part and n is the absolute value 
of the imaginary part of the remaining two complex conjugate roots of the polynomial in Eq.(^5|). 

In the weak field case there are 4 special values of the Hamiltonian, Hmin, the minimum value, H saddle the value 
at the saddle point, Hmaxi and Hmax2, the smaller and the larger of the values at the two maxima respectively (See 
Figure 1). For H values in the range Hmax2 > H' > Hmaxi or H saddle > H > Hmin there is again one to one 
correspondence between H values and trajectories on the sphere. For such trajectories the frequencies are again given 
by (|4^). For the values of H satisfying Hmaxi > H > H saddle there are two trajectories corresponding to each H, one 
circling around Hmaxi, and the other circling around Hmax2- However they both have the same period given by 



T= 4 {X2-X,)[x,^x,) ^^^^ 

\Ay{xi - X2){X3 - Xi) y [Xi - X2)[X3 - Xi) 

where K is again the complete elliptic integral and xi < X2 < x^ < Xi are the four real roots of the polynomial in 
Eq.®. 

A typical plot of frequencies for the weak field case is given in figure 2. Near the saddle point trajectories slow 
down logarithmically in _ff — H saddle as expected from a two dimensional dynamics. 
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We can see for An = A22 that our equations reduce to those obtained in |3j:^l[l by using two coupled Gross- 
Pitaevskii equations. At high field we get almost circular trajectories, and correspondingly large oscillations between 
internal levels of the condensate, as in the Rabi problem. For the low field case we get three kind of trajectories all 
of which give little population change, these correspond to Josephson oscillations. Two of these three kinds complete 
a full cycle around the z axis, while the third is trapped in a region for which (pmin < < (t>max- Recalling that (j) 
represents the relative phase of the two condensates, we see that these trajectories correspond to Josephson oscillations 
between the two internal states, caused by the weak link of the external field. The other two kinds of trajectories 
again correspond to Josephson oscillations, however in these class of Josephson oscillations there is a 2tt phase slip 
for every period of population change, such oscillations are called tt Josephson oscillations [Q. 

It has been shown for the spatial Josephson effect that the system can dynamically maintain a population imbalance 
between the two macroscopically occupied states ^2j. We can ask whether the internal Josephson effect also shows 
this "macroscopic quantum self trapping" effect. This question is meaningful only if a symmetry exists between states 
1 and 2 

Aii = A22. (49) 

If the internal states have different interactions a population imbalance is to be expected in any oscillation as the 
system will be biased to choose the internal state with lower interaction energy. 

If we take An = A22, our condition for the appearance of the extra maximum in the Hamiltonian becomes 

|C|>2|^|, (50) 

as -B = in this case. One can immediately see that out of the three types of trajectories discussed for the general 
case, those corresponding to tt Josephson oscillations become macroscopically self trapped oscillations for An = A22. 

It is remarkable that none of the arguments given for the internal Josephson effect depend on the macroscopic 
occupation of a spatial state, i.e. Bose condensation. In fact if we assume that the system is not condensed 

T > Te, (51) 

we get the equation for the precession of the normal spin from Eq.([21 



S^Sx Bn + 2SxV (52) 

Bn = [(An - A22)2/9o + (An + A22 - 2Ai2)25" • z]z. 

This equation is very similar to the equation for the condensate spin, and all the results found in this section apply 
for the non-condensed case if we identify 

^ = (An+A22-2Ai2)|5p (53) 
B = i2{Xn-X22)po\S\ + 2V'\S\) 
C = 2^131 
D = 2X,2\S\\ 

The behavior of the fully condensed and non-condensed cases are similar, however the values of the precession fre- 
quencies and critical value of the external field are different. This difference, again, is a result of the fact that there 
are no exchange scattering processes in a Bose condensate. 

From the above discussion we come to the conclusion that to observe the internal Josephson effect it is not required 
to have two Bose condensed samples. Two words of caution should be voiced about this observation. First we have 
only considered the coherent collisions and the coherence of the phase of two internal levels will be destroyed on a 
time scale that is set by mean free time in the gas. Any observed internal Josephson oscillation in a non-condensed 
sample should decay in this time scale. Second, to be able to see this effect one has to go to very low field strengths 
such that the population oscillations should be observable in a non-condensed sample. 

Still we have shown that it is not absolutely necessary to have Bose condensation to observe small oscillations in 
the relative phase of the internal states under a small excitation field. Let us imagine a non-condensed gas of atoms 
put into a superposition of two internal states by a 7r/2 pulse. If this sample is further subjected to weak mixing 
field on resonance one would naively expect population transfer from one internal level to the other with the Rabi 
frequency of the field. However our discussion shows that if the sample and the field satisfy 

T^oi\ = SnalnvT < \V\ < An (54) 

we would have small oscillations of the phase and populations, which is exactly what is observed in the internal 
Josephson effect with condensed samples. 
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VI. ABSORPTION RESONANCES FOR A PARTIALLY CONDENSED GAS: THE TWO SPIN PROBLEM 



After analyzing the fully condensed and non-condensed Bose gas, both of which can be represented by only one spin 
in our equations (|2^), we turn our attention to the partially condensed Bose gas. When the temperature is between 
zero and condensation temperature, TbeC) there is both a condensed and a non-condensed density present. If these 
densities are comparable we have to use the full form of Eq.(^7|), with both spins present. 

We first start with the free precession problem, i.e. by setting the external field equal to zero. In this case we have 
the equations 

S = Sx B„ + 2X12S X Sc (55) 
Sc^ScX Bc + 2\i2Sc X S. (56) 



The effective magnetic fields Bn and Be both point in the z direction and are defined in Eq.(^8|). As their time deriva- 
tives are perpendicular to the spin vectors, norm of both of the spins, \S\ and \Sc\ are conserved. This conservation 
simply means that the total number of atoms in the condensate and over the condensate are conserved. Another 
conserved quantity can be obtained by adding Eq.(|55|) to Eq.(|5^) and taking the z component. We have 

|(i.(5 + 5.)).|54, = 0. (57) 

Physically this conservation law corresponds to the fact that in the absence of a coupling field the total number 
of atoms in internal states 1 and 2 are conserved separately. Although these three conservation laws restrict the 
resulting dynamics considerably, they still allow oscillations in which the density of condensed and non-condensed 
atoms in internal state 1 change, while state two goes through the same oscillations out of phase with 1 to leave the 
total number of condensed and non-condensed atoms constant. In our spin representation the degree of freedom that 
expresses these oscillations will be S*^, or equivalently 5f as they add up to a constant. 

The other conserved quantity is the Hamiltonian ( ^ ) which we choose to rewrite in terms of the conserved quantities 
S't^ot , \S\ , \Sc\ and the dynamical variable S^ as 

n = (An - A22)(po + Po)S^ot + ^(Aii + A22 - 2X,2)iSU' (58) 
+ 2Ai2|5|2 + Ai2|5,p + (An - A22)po5" 
+ ^(An + A22 - 2Ai2)(5^)2 + 2Ai25 • S^. 

To investigate the oscillations of the degree of freedom physically described above and represented by S^ we take 
the equation of motion for S^ 

— = z-S = 2A122 • (S X Sc) - 2A12V. (59) 
dt 

We define V to be the volume of the parallelepiped formed by the vectors z , S and Sc. We can express the absolute 
value of this volume in terms of the inner products of these three vectors as 

|V| = v^|5|2|5c|2 - (S-Scy - |5e|2(5-)2 - \S\^ {S + 2S^ S^ S ■ Sc. (60) 

Now we can solve for S* • 5c in terms of conserved quantities and 5^ from Eq.(p8[). This will give us the equation of 
motion for 5^ expressed only in terms of conserved quantities and 5^ itself: 



dS' 



dt 



= VC4(5-)4 + CsiS^r + C2(5-)2 + Ci{S-) + Co (61) 



The value of 5^ will oscillate between xi and X2 which are two roots of the polynomial inside the square root in 
Eq. (IsTl) . For values of 5^ in the interval xi < S^ < X2, this polynomial takes positive values. When 5^ reaches its 
maximum or minimum value, the vectors 5 , Sc and z are coplanar. This allows us to integrate Eq. (|6l|) without 
paying attention to the absolute value. We can express the period of 5^ as an integral of the form in Eq.(|45|) 

l'X2 I 

Tz = 2 dS^ ,^ ^ ^ ^ _ ^ . (62) 



VC4(5^)4 + C3(5^)3 + C2(5^)2 + CiS^ + Co ■ 
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Due to the abundance of conserved quantities in the two spin problem, the expressions for the coefficients are more 
complicated compared to the one spin case 



C3 
C2 

Ci 
Co 



A2 A2 
2 ^ 2 



4Ai 



(63) 



[AiA2Po + 4Ai2(^^tot-AiPo)] 



2Ai2|^P + Ai2|5,n 
4Ai2(H - Aipo5f„t - ^{S^,,r + Ai2|^n] 
{■H ~ Ai(po + Po)S^ot - ^{SLt? - 2Ai2|5p + Ai2|5cn(2Aipo + ^\i2Slt) 



%x\2\s\^si 



- (n - Ai(po + Po)S^ot - -Yi^totf - 2Ai2|5p + Ai2|^cn' 
+ 2in - Ai(po + po)^t^t)(2Ai2|5p + Ai2|5,n ~ 4A?2|^|* 
- AL|5.|4 - A2Xi2\Scf{S^,,f - 2(A2 + 2Ai2)Ai2|5|2(54t)', 

with the notation 

Ai = (All — A22) 

A2 = (All + A22 — 2A12). 



(64) 



As in the one spin case this integral can be exactly evaluated 
square root in Eq.(|6l|) are real we have 



If all the roots of the polynomial inside the 



C4 y/{x4 - X2){X3 - Xl) 



-.K 



[X2 - Xl){x4 - X3) 

(0:4 - X2){x:i - Xl) 



(65) 



Here X3 < X4 are the remaining real root of the polynomial in equation of motion (^Tj), which are assumed to be real. 
In the case of imaginary 2:3 and X4, an analogue of Eq.(^) will give the expression for the period. 

The second oscillation for free precession corresponds to the precession of total phase about the z axis. This pre- 
cession is affected by the population oscillations found above, and the oscillation frequency is not as easily calculated. 
We can generally describe its motion as a sum of two components. One corresponding to a uniform precession around 
the z axis with the density shift as in Eq(^) , and the other corresponding to the effect of an oscillating magnetic field 
in z direction, caused by the population oscillations discussed above. The coupling between these two components is 
best seen when we write the equation of motion for the spin components in the x — y plane. Defining = S*^ + iS"^ 
and S+ = 8"^ + iSy we have 



iS+ = (B„(t) + 2Xi2S'^{t))S+ - 2Xi2S'{t)S+ 
iS+ = -2Ai25,"(t)5+ + {Be + 2Xi2S'it))S+. 



(66) 



We have seen that when reaches its maximum or minimum values S and Sc are in the same vertical plane. We 
can calculate exactly how much the spins have rotated around the z axis throughout the course of one oscillation. 
The effect of oscillations will present itself through the integral 



dtS'{t) 



X2 QZ 



^1 I dt 



(67) 



which is again exactly calculable. 

To calculate the rotation angle we integrate Eq. (p&j . After time Tz , 5*+ and 5*+ will be given by 







■ ^+(0) " 









(68) 



where is a two by two matrix with elements 
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Mn = Ai(2po + Po)T, + (A2 + 2Xi2)S^^,T, + (A2 - 2Ai2)/ (69) 

M12 = -2A12/ 

M21 = -2\i2Sl^T, + 2A12/ 

M22 = Ai(po + Po)T, + ^2Sl,T, + 2A12/. 

Instead of giving the resulting long expression for the rotation angle in one period, we choose to describe the motion 
qualitatively. The precession of the two spins are affected by the competition between two effects. The first effect is, 
due to the absence of exchange scattering in the condensate, the condensate spin Sc sees an effective magnetic field 
Be different from the effective magnetic field seen by the normal gas spin S. The second, as discussed above, is the 
condensate population oscillations characterized by 5^. 

If there is not much difference between the densities of two internal states, both spins lie close to the x-y plane, 
and their relative phase oscillates around zero, never growing large. However, if there is a lot of density difference 
between two internal states, the spins are close to the z axis. Over one period of S*^ oscillation, the phase difference 
can be a multiple of 27r. 

We can investigate the precession easily in the limit when both of the spins are almost aligned with the z axis. We 
can write linear equations for the perpendicular components of the spins and get the two precession frequencies: 

w = a) + (w2 - wi) + (Ai2(ni - 712) + 2A22n2 - 2Aiini). (70) 

Where a) satisfies, 

[Co + Ai2A7it)(a) - AiiTici + A22nc2) = Ai2AncAnt. (71) 

with 

An(c,t) = n2(c,t) - ni(c^t). (72) 

If the mixing angle is not small one would expect to see two different frequencies in general. The average of the 
two frequencies will be controlled by the average density shift seen by an atom in the sample, while the splitting will 
reflect the average rate the condensate fraction of one of the internal states oscillates p3| . 

The appearance of a second frequency should be detectable in an experiment that probes a partially condensed Bose 
gas in a superposition state. We propose using a partially condensed gas in a Ramsey separated field arrangement, 
as in the fountain atomic clocks [ p^|jl7| ]. In this case the appearance of a second frequency would present itself as a 
beating in the Ramsey fringes. This beating however, will vanish both in the limit of full condensation and the limit 
of a normal gas, and should be most prominent when the condensate fraction is close to a half. The exact values of 
the frequencies can be obtained solving the equations of motion numerically. 

Here we want to remind the reader that the equations used in this section Eq. (|2^) were derived for a uniform system. 
If the particles are cold enough and the condensate is prepared in a shallow trap to make sure that the movement 
of each part of the cloud is negligible in the center of mass coordinate frame during the time of measurement, the 
equations will be locally correct and the experiment should show a density averaged result in the precession frequencies. 
Otherwise the effects of inhomogeneity must be included using the full transport equations Eq.(17 T^). 



When an external field is turned on the z component of total spin in not conserved anymore, there is a net population 
transfer from one internal state to the other. Much like the single spin case we have two limits. When \V\ is much 
larger then the density shift An we have both spins following almost circular trajectories around V . For the weak 
field case the system can be best described as two non-linear oscillators going through coupled oscillations with an 
occasional two tt phase slip for one of them. In this most general case precession frequencies can be found by numerical 
integration of Eq.(p7|). 



VII. FERMIONS 



Finally to understand the effect of statistics better we consider the same problem for a Fermi gas. We will have the 
same Hamiltonian Eq.(^, however, the Fermionic field operators will satisfy 

{Mr),i>+p{r')}^5a.p5{r~r'), (73) 

where {,} denotes the anti-commutator. The derivation will process along the same lines with the Bose case. The 
effect of statistics will be seen whenever we average a four particle operator. The exchange term in Eq.(|l6|) will change 
sign: 
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(74) 



As a result we get the transport equation for the density matrix defined as in Eq.(^^ with sign changes in the 
terms corresponding to the exchange contributions: 



ft + |-v.-v.(Mi±Ml).v,),„,(.,p) 



(75) 



Ua{r) ~ Up{r) + ^(A-ya - A-y^)p^^(r) 



E 



(r) • Vpg^,3{r,p) + ^ Xp^VrP^pir) ■ {r,p) 

7 7 

+ *E 1^7" t)Q-(p{r,p) - Vp^{r, t)gay{r,p)] 

7 

{r,p)]. 

7 

Prom the transport equation by assuming all the interactions and the sample to be spatially homogeneous we can 
get the equation of motion for the internal state density matrix 



with are defined as, 



^77' ^ ('^^7 ' I J-J'^ t^-y ' ) /^"yY' ^ ^ (A(^-y X(^^^^ P^(^ Pq 

a 

+i^(VaT,(t)pa7' " VYa{t)pja) 



Ua = ^ Xapnj3. 



(76) 



(77) 



If we assume that only states 1 and 2 are coupled and all the off diagonal elements involving the other states are 
equal to zero, we get a very simple dynamics. The time derivative of the diagonal elements do not depend on p while 
the off diagonal element pi2 changes according to 



P12 = l \ Wi- W2+ ^ (Xf3l - Xj32)np j P12 



(78) 



+i 



'^(Val{t)pa2 - V2a{t)pia) 



When we go to the Bloch sphere representation in the basis rotating with the frequency of the external field: 

p'-yY = PoS^Y + Sf ■ a^Y C^S) 

we get the simple equation of motion 

Sf = 2SfxV. (80) 
Which can be derived from the corresponding Hamiltonian 

Hf = 2Sf ■ V. (81) 
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In the mean field picture there are no coherent effects of interactions for a transition between two states. The 
exchange contributions to the precession frequency exactly cancel the direct contributions. For short-range potentials, 
the exchange contribution to the energy has the same absolute value as the direct contribution. For Bosons the two 
contributions add up, while for fermions they cancel. 

The precession frequency for a free Fermi gas can be read off from Eq. ( fz^ ) 

u = [w2 - wi) + ^ (A2t- - \i^)n~f. (82) 

This expression shows that the density dependent frequency shift encountered in the fountain atomic clocks can be 
eliminated if a fermionic sample is used instead of a Bose gas. Any contributions to the frequency shift will be of 
higher order in the dilutcness parameter as/{n~^/^) of the gas. Finally we see that the behavior under an external 
filed is not at all different from the usual Rabi precession. An analogue of internal Josephson effect does not appear 
in this case since the "energy" of the system Eq.(^) does not depend on the density at all (there are no quadratic 
terms in the Hamiltonian in the spin representation). 



VIII. CONCLUSION 



We studied the effect of external interactions on the internal dynamics of atoms in a dilute gas. For a Bose gas 
we derived a general transport equation valid for any partially condensed and/or non-uniform gas. We then applied 
it to the case of a homogeneous gas and investigated the effects of interactions on the internal degrees of freedom. 
As a first result we obtained an expression for the density induced frequency shift in atomic clocks, for a gas which 
is above BEC, or is at zero temperature. Furthermore we found that if a partially condensed sample is used in an 
atomic clock, one would get two density dependent frequencies instead of one, due to the exchange of atoms between 
the normal part and the condensed part of the gas. 

We then went on to analyze the effect of an external field. We show how Rabi oscillations are replaced by internal 
Josephson oscillations as the strength of the external field is reduced. We calculated the frequencies of both oscillations 
exactly. We have also found that an analogue of the internal Josephson effect should be observable for a non-condensed 
sample. 

Finally we considered a Fermi gas, and derived the transport equation. We demonstrated that it is possible to 
eliminate density shift in Rabi frequency by using a Fermi gas in an atomic clock and that no analogue of the internal 
Josephson effect is possible for a Fermi gas. 

M.0.0. is grateful to S.R. Shenoy, F. Sols and D. Stroud for useful discussions. 

t Current address: Department of Physics, The Ohio State University, Columbus OH 43210. 
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FIG. 1. The contour lines for the Hamiltonian Eq.(^0|), on the Bloch sphere, defined by , 0, for weak external fields. 
As the field strength is increased, the saddle point comes closer to one of the maxima and they destroy each other when the 
condition in Eq.(k3) is satisfied, leaving just one maximum and one minimum. In the figure the trajectories encircling the 
maxima correspond to Josephson oscillations with an average phase difference of n, while those encircling the minimum are 
the usual Josephson oscillations corresponding to small oscillations of phase difference. The frequencies of motion along these 
trajectories are given in fig2. 




H (Arbitrary Units) l_. 
FIG. 2. The frequency of precession on the Bloch sphere as a function of the value the Hamiltonian Eq.(WO|) takes. 

Hamiltonian can take values from Hmin to Hmax2- Near the saddle point the precession slows down logarithmically, as 

expected from a two dimensional dynamics. Between Hsaddie and Hmaxi there are two trajectories for each value that the 

Hamiltonian takes. However, they both have the same frequency Eq.([l8|). 
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